Search results

The main purpose of this paper is to compare the performance of three commonly used global search algorithms, namely tree‐based augmented spatial digital tree, cell‐based no binary search and D‐cell, in the discrete element simulations.

A large number of test cases with up to five million particles/discrete objects are employed to numerically examine the computational costs of the three search algorithms and their performance is compared.

Comprehensive comparisons reveal that the D‐cell is more efficient than the tree‐based search algorithms for large‐scale problems. The parametric study of the D‐cell algorithm itself shows that the performance of the algorithm is strongly dependent on the cell dimension chosen.

The only limitation of the current work is that the tested domain shape is regular, and thus more complex domain shapes may need to be considered.

The paper provides clear guidance regarding the possible actual computational performance of the tested search algorithms for practical applications.

3D steady heat conduction analysis considering heat source is conducted on the fundamental of the fast multipole method (FMM)-accelerated line integration boundary element method (LIBEM).

Due to considering the heat source, domain integral is generated in the traditional heat conduction boundary integral equation (BIE), which will counteract the well-known merit of the BEM, namely, boundary-only discretization. To avoid volume discretization, the enhanced BEM, the LIBEM with dimension reduction property is introduced to transfer the domain integral into line integrals. Besides, owing to the unsatisfactory performance of the LIBEM when it comes to large-scale structures requiring massive computation, the FMM-accelerated LIBEM (FM-LIBEM) is proposed to improve the computation efficiency further.

Assuming N and M are the numbers of nodes and integral lines, respectively, the FM-LIBEM can reduce the time complexity from O(NM) to about O(N+ M), and a full discussion and verification of the advantage are done based on numerical examples under heat conduction.

(1) The LIBEM is applied to 3D heat conduction analysis with heat source. (2) The domain integrals can be transformed into boundary integrals with straight line integrals by the LIM. (3) A FM-LIBEM is proposed and can reduce the time complexity from O(NM) to O(N+ M). (4) The FM-LIBEM with high computational efficiency is exerted to solve 3D heat conduction analysis with heat source in massive computation successfully.

The purpose of this paper is to extend the previous study [Computer Methods in Applied Mechanics and Engineering 340: 70-89, 2018] on the development of a novel packing characterising system based on principal component analysis (PCA) to quantitatively reveal some fundamental features of spherical particle packings in three-dimensional.

Gaussian quadrature is adopted to obtain the volume matrix representation of a particle packing. Then, the digitalised image of the packing is obtained by converting cross-sectional images along one direction to column vectors of the packing image. Both a principal variance (PV) function and a dissimilarity coefficient (DC) are proposed to characterise differences between different packings (or images).

Differences between two packings with different packing features can be revealed by the PVs and DC. Furthermore, the values of PV and DC can indicate different levels of effects on packing caused by configuration randomness, particle distribution, packing density and particle size distribution. The uniformity and isotropy of a packing can also be investigated by this PCA based approach.

Develop an alternative novel approach to quantitatively characterise sphere packings, particularly their differences.

Contact detection for convex polygons/polyhedra has been a critical issue in discrete/discontinuous modelling, such as the discrete element method (DEM) and the discontinuous deformation analysis (DDA). The recently developed 3D contact theory for polyhedra in DDA depends on the so-called entrance block of two polyhedra and reduces the contact to evaluate the distance between the reference point to the corresponding entrance block, but effective implementation is still lacking.

In this paper, the equivalence of the entrance block and the Minkowski difference of two polyhedra is emphasised and two well-known Minkowski difference-based contact detection and overlap computation algorithms, GJK and expanding polytope algorithm (EPA), are chosen as the possible numerical approaches to the 3D contact theory for DDA, and also as alternatives for computing polyhedral contact features in DEM. The key algorithmic issues are outlined and their important features are highlighted.

Numerical examples indicate that the average number of updates required in GJK for polyhedral contact is around 6, and only 1 or 2 iterations are needed in EPA to find the overlap and all the relevant contact features when the overlap between polyhedra is small.

The equivalence of the entrance block in DDA and the Minkowski difference of two polyhedra is emphasised; GJK- and EPA-based contact algorithms are applied to convex polyhedra in DEM; energy conservation is guaranteed for the contact theory used; and numerical results demonstrate the effectiveness of the proposed methodologies.

The purpose of this paper is to establish a new two-phase Discrete Element Method (DEM) model to investigate the movement of fresh concrete which consists of mortar and aggregate. The established DEM model was adopted to simulate the mixing process of fresh concrete based on the commercial software package PFC3D. The trajectories of particles and particle clusters were recorded to analyze the mixing behavior from different scales. On one hand, the macro-scale movement was obtained to make the mixing process visualization. On the other hand, the relative micro movement of the single particle and particle clusters was also monitored to further study the mixing mechanism of the fresh concrete.

A new two-phase DEM model was designed to simulate the movement of fresh concrete which consists of mortar and aggregate. The linear-spring dashpot model was used to model all the contacts between particle and particle/wall to characterize the viscidity of fresh concrete. Moreover, two sets of parallel bond models were employed to characterize the contact between the mortar particles and mortar/coarse aggregate particles, namely the pbond1 and pbond2. The hybrid treatment enables the current DEM model to handle the yield behavior.

The mixing process of fresh concrete is mainly composed by the transportation in the x-direction and the overturn and fall off in the y- and z-directions. With these movements in different directions, the concrete particles can be fully mixed in the mixing drum.

A new two-phase DEM model was proposed and used to simulate the mixing process of fresh concrete. The outcomes of the simulation would be helpful for making the transporting truck visualization and the movement behavior of fresh concrete observable. The model can provide dynamic information of particles to reveal the interaction mechanism of fresh concrete in the truck mixer which is extremely difficult to obtain on-line in physical experiments or building site.

The purpose of this paper is to improve the comfort and shaping function of seamless shapewear on the material and structure and develop new seamless shapewear products. The shaping figure effect will be verified as well.

The performance of the knitted fabrics made of Polytrimethylene terephthalate (PTT) filament was analyzed by orthogonal experiment and fuzzy mathematical methods analysis, in order to get the optimal conditions for the best performance. The new products were designed and made based on the results of the material research with the consideration of the aesthetic requirements. The shaping effect of seamless shapewears on local and global figure was tested by the methods of the combination of subjective and objective evaluation.

The sample which renders the optimal performance for shapewear is the one with PTT filament as face yarn, nylon core-spun yarn as ground yarn and 3+1 simulate rib knit structure. The material of face yarn, mixed proportion and structure can influence the shape retention, appearance and comfort of PTT fabric in various degrees. Three shapewears which were developed according to the results of material research have different shaping effect. And women with different figures put different satisfaction degrees on each shapewear’s shaping effect.

This paper provides scientific basis and reference for enterprise to design good tight seamless shapewear as well as for consumer to buy suitable products.

In the view of the problems of present shapewears, this paper completed the development of the shapewears and verified the shaping effect of them on women with different figures. The shapewears can be put into production directly.

This study uses migrant household survey data from 2008 to 2009 to examine how parental migration decisions are associated with the nutritional status of children in rural and urban China. Results from instrumental variables regressions show a substantial adverse effect of children’s exposure to parental migration on height-for-age Z scores of left-behind children relative to children who migrate with their parents. Additional results from a standard Blinder–Oaxaca decomposition, a quantile decomposition, and a counterfactual distribution analysis all confirm that children who are left behind in rural villages – usually because of the oppressive hukou system – have poorer nutritional status than children who migrate with their parents, and the gaps are biggest at lower portions of the distribution.

Construction insurance has been advocated as a useful financial means to transfer risks and mitigate potential project losses. However, there is a general reluctance to purchase construction insurance in China. The purpose of this paper is to develop an extended theory of planned behavior (TPB) model to understand and predict contractors’ construction insurance purchasing intention.

Quantitative data were collected through a questionnaire survey. The respondents comprise 130 project directors/managers and contract managers/administrators who were involved in or familiar with the construction insurance purchasing or managing. Data were analyzed using structural equation modeling.

This study developed an extended TPB model to explain the contractors’ intention formation. The results indicated that the insurance purchasing intention was mainly influenced by attitudes and subjective norm, rather than perceived behavioral control. The analyses also revealed that the added variables (risk perception and past experience) had a significant impact on the attitudes and perceived behavioral control.

This study identified the factors leading to Chinese contractors’ low willingness to purchase construction insurance by extending the TPB model to the context of construction insurance purchasing. The extended TPB model may serve as a theoretical framework and basis for the Chinese contractors, insurers, and regulators to understand the root causes of contractors’ low willingness to purchase construction insurance and make joint efforts to address the risks and insurance in the construction industry.

Sand production is a challenging issue during hydrocarbon production in the oil and gas industry. This paper aims to investigate one sand production process, i.e. transient sand production, using a novel bonded particle lattice Boltzmann method. This mesoscopic technique provides a unique insight into complicated sand erosion process during oil exploitation.

The mesoscopic fluid-particle coupling is directly approached by the immersed moving boundary method in the framework of lattice Boltzmann method. Bonded particle method is used for resolving the deformation of solid. The onset of grain erosion of rocks, which are modelled by a bonded particle model, is realised by breaking the bonds simulating cementation when the tension or tangential force exceeds critical values.

It is proved that the complex fluid–solid interaction occurring at the pore/grain level can be well captured by the immersed moving boundary scheme in the framework of the lattice Boltzmann method. It is found that when the drawdown happens at the wellbore cavity, the tensile failure area appears at the edge of the cavity. Then, the tensile failure area gradually propagates inward, and the solid particles at the tensile failure area become fluidised because of large drag forces. Subsequently, some eroded particles are washed out. This numerical investigation is demonstrated through comparison with the experimental results. In addition, through breaking the cementation, which is simulated by bond models, between bonded particles, the transient particle erosion process is successfully captured.

A novel bonded particle lattice Boltzmann method is used to investigate the sand production problem at the grain level. It is proved that the complex fluid–solid interaction occurring at the pore/grain level can be well captured by the immersed moving boundary scheme in the framework of the lattice Boltzmann method. Through breaking the cementation, which is simulated by bond models, between bonded particles, the transient particle erosion process is successfully captured.